Evaluation of oyster exposure experiments
knitr::opts_chunk$set(echo = TRUE, warning = F, message = F, fig.path = 'figs/', dev.args = list(family = 'serif'))
library(tidyverse)
library(patchwork)
library(ggbeeswarm)
library(psycho)
library(lmerTest)
library(hrbrthemes)
library(kableExtra)
library(shiny)
source('R/funcs.R')
data(allexp)
data(bymods)
data(wtmods)
Number of folds
Raw data
yvar <- 'Number of folds'
toplo <- allexp %>%
filter(var %in% yvar)
p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) +
geom_boxplot(outlier.shape = NA, alpha = 0.1) +
geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) +
theme_ipsum_rc() +
facet_grid(trt~ .) +
scale_colour_ipsum() +
scale_fill_ipsum() +
theme(legend.title = element_blank()) +
labs(x = 'Exposure week', y = yvar)
p
Models with week
Olympia oyster
dat <- wtmods %>%
filter(species == 'Olympia') %>%
filter(var == yvar)
dat %>%
filter(species == 'Olympia') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Number of folds vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.74 |
0.18 |
4 |
32.45 |
11.48 |
0.00 |
| week |
0.01 |
0.01 |
2 |
30.48 |
0.38 |
0.69 |
dat %>%
filter(species == 'Olympia') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Number of folds vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.63 |
0.16 |
4 |
25.69 |
9.92 |
0.00 |
| week |
0.01 |
0.01 |
2 |
24.16 |
0.42 |
0.66 |
| trt:week |
0.07 |
0.01 |
8 |
25.46 |
0.58 |
0.79 |
dat %>%
pull(plomod)
Pacific oyster
dat <- wtmods %>%
filter(species == 'Pacific') %>%
filter(var == yvar)
dat %>%
filter(species == 'Pacific') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Number of folds vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.66 |
0.17 |
4 |
38.14 |
8.31 |
0.00 |
| week |
0.01 |
0.00 |
2 |
36.80 |
0.21 |
0.81 |
dat %>%
filter(species == 'Pacific') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Number of folds vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.73 |
0.18 |
4 |
31.25 |
9.19 |
0.00 |
| week |
0.01 |
0.01 |
2 |
30.19 |
0.36 |
0.70 |
| trt:week |
0.23 |
0.03 |
8 |
31.05 |
1.42 |
0.23 |
dat %>%
pull(plomod)
Number of irregularites
Raw data
yvar <- 'Number of irregularities'
toplo <- allexp %>%
filter(var %in% yvar)
p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) +
geom_boxplot(outlier.shape = NA, alpha = 0.1) +
geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) +
theme_ipsum_rc() +
facet_grid(trt~ .) +
scale_colour_ipsum() +
scale_fill_ipsum() +
theme(legend.title = element_blank()) +
labs(x = 'Exposure week', y = yvar)
p
Models with week
Olympia oyster
dat <- wtmods %>%
filter(species == 'Olympia') %>%
filter(var == yvar)
dat %>%
filter(species == 'Olympia') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Number of irregularities vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.31 |
0.08 |
4 |
37.07 |
1.68 |
0.18 |
| week |
0.06 |
0.03 |
2 |
31.63 |
0.64 |
0.54 |
dat %>%
filter(species == 'Olympia') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Number of irregularities vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.38 |
0.10 |
4 |
26.01 |
2.24 |
0.09 |
| week |
0.07 |
0.03 |
2 |
23.62 |
0.80 |
0.46 |
| trt:week |
0.53 |
0.07 |
8 |
26.22 |
1.54 |
0.19 |
dat %>%
pull(plomod)
Pacific oyster
dat <- wtmods %>%
filter(species == 'Pacific') %>%
filter(var == yvar)
dat %>%
filter(species == 'Pacific') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Number of irregularities vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.45 |
0.11 |
4 |
38.34 |
4.00 |
0.01 |
| week |
0.02 |
0.01 |
2 |
36.72 |
0.32 |
0.73 |
dat %>%
filter(species == 'Pacific') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Number of irregularities vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.70 |
0.18 |
4 |
27.51 |
6.08 |
0.00 |
| week |
0.03 |
0.02 |
2 |
26.70 |
0.53 |
0.60 |
| trt:week |
0.78 |
0.10 |
8 |
27.19 |
3.40 |
0.01 |
dat %>%
pull(plomod)
Total observations (dissolution)
Raw data
yvar <- 'Total observations'
toplo <- allexp %>%
filter(var %in% yvar)
p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) +
geom_boxplot(outlier.shape = NA, alpha = 0.1) +
geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) +
theme_ipsum_rc() +
facet_grid(trt~ .) +
scale_colour_ipsum() +
scale_fill_ipsum() +
theme(legend.title = element_blank()) +
labs(x = 'Exposure week', y = yvar)
p
Models with week
Olympia oyster
dat <- wtmods %>%
filter(species == 'Olympia') %>%
filter(var == yvar)
dat %>%
filter(species == 'Olympia') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Total observations vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
1.16 |
0.29 |
4 |
39.12 |
1.19 |
0.33 |
| week |
0.14 |
0.07 |
2 |
35.48 |
0.30 |
0.75 |
dat %>%
filter(species == 'Olympia') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Total observations vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
1.42 |
0.36 |
4 |
31.93 |
1.48 |
0.23 |
| week |
0.16 |
0.08 |
2 |
29.77 |
0.33 |
0.72 |
| trt:week |
2.37 |
0.30 |
8 |
31.83 |
1.23 |
0.31 |
dat %>%
pull(plomod)
Pacific oyster
dat <- wtmods %>%
filter(species == 'Pacific') %>%
filter(var == yvar)
dat %>%
filter(species == 'Pacific') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Total observations vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
1.55 |
0.39 |
4 |
36.54 |
1.16 |
0.34 |
| week |
2.40 |
1.20 |
2 |
34.90 |
3.61 |
0.04 |
dat %>%
filter(species == 'Pacific') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Total observations vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
1.87 |
0.47 |
4 |
28.70 |
1.40 |
0.26 |
| week |
2.71 |
1.36 |
2 |
27.63 |
4.07 |
0.03 |
| trt:week |
4.51 |
0.56 |
8 |
28.56 |
1.69 |
0.14 |
dat %>%
pull(plomod)
Dissolution score
Raw data
yvar <- 'Dissolution score (0-3)'
toplo <- allexp %>%
filter(var %in% yvar)
p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) +
geom_boxplot(outlier.shape = NA, alpha = 0.1) +
geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) +
theme_ipsum_rc() +
facet_grid(trt~ .) +
scale_colour_ipsum() +
scale_fill_ipsum() +
theme(legend.title = element_blank()) +
labs(x = 'Exposure week', y = yvar)
p
Models with week
Olympia oyster
dat <- wtmods %>%
filter(species == 'Olympia') %>%
filter(var == yvar)
dat %>%
filter(species == 'Olympia') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Dissolution score (0-3) vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
3.90 |
0.97 |
4 |
39.04 |
9.20 |
0.00 |
| week |
0.07 |
0.03 |
2 |
36.14 |
0.32 |
0.73 |
dat %>%
filter(species == 'Olympia') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Dissolution score (0-3) vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
4.78 |
1.19 |
4 |
28.14 |
10.96 |
0.00 |
| week |
0.05 |
0.02 |
2 |
26.69 |
0.22 |
0.81 |
| trt:week |
1.37 |
0.17 |
8 |
27.81 |
1.58 |
0.18 |
dat %>%
pull(plomod)
Pacific oyster
dat <- wtmods %>%
filter(species == 'Pacific') %>%
filter(var == yvar)
dat %>%
filter(species == 'Pacific') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Dissolution score (0-3) vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
3.24 |
0.81 |
4 |
37.38 |
10.71 |
0.0 |
| week |
0.14 |
0.07 |
2 |
36.11 |
0.95 |
0.4 |
dat %>%
filter(species == 'Pacific') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Dissolution score (0-3) vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
3.30 |
0.82 |
4 |
27.73 |
10.76 |
0.00 |
| week |
0.15 |
0.07 |
2 |
26.72 |
0.97 |
0.39 |
| trt:week |
0.49 |
0.06 |
8 |
27.56 |
0.80 |
0.61 |
dat %>%
pull(plomod)
Shell weight
Raw data
yvar <- 'Shell weight (g)'
toplo <- allexp %>%
filter(var %in% yvar)
p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) +
geom_boxplot(outlier.shape = NA, alpha = 0.1) +
geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) +
theme_ipsum_rc() +
facet_grid(trt~ .) +
scale_colour_ipsum() +
scale_fill_ipsum() +
theme(legend.title = element_blank()) +
labs(x = 'Exposure week', y = yvar)
p
Models with week
Olympia oyster
dat <- wtmods %>%
filter(species == 'Olympia') %>%
filter(var == yvar)
dat %>%
filter(species == 'Olympia') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Shell weight (g) vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.04 |
0.01 |
4 |
61.89 |
0.30 |
0.88 |
| week |
0.10 |
0.05 |
2 |
60.58 |
1.36 |
0.26 |
dat %>%
filter(species == 'Olympia') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Shell weight (g) vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.05 |
0.01 |
4 |
51.80 |
0.32 |
0.86 |
| week |
0.10 |
0.05 |
2 |
52.45 |
1.29 |
0.28 |
| trt:week |
0.26 |
0.03 |
8 |
53.41 |
0.86 |
0.56 |
dat %>%
pull(plomod)
Pacific oyster
dat <- wtmods %>%
filter(species == 'Pacific') %>%
filter(var == yvar)
dat %>%
filter(species == 'Pacific') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Shell weight (g) vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.05 |
0.01 |
4 |
55.71 |
1.29 |
0.29 |
| week |
0.03 |
0.01 |
2 |
53.03 |
1.29 |
0.28 |
dat %>%
filter(species == 'Pacific') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Shell weight (g) vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.06 |
0.01 |
4 |
44.24 |
1.45 |
0.23 |
| week |
0.03 |
0.01 |
2 |
43.73 |
1.35 |
0.27 |
| trt:week |
0.11 |
0.01 |
8 |
44.50 |
1.41 |
0.22 |
dat %>%
pull(plomod)
Tissue weight
Raw data
yvar <- 'Tissue weight (g)'
toplo <- allexp %>%
filter(var %in% yvar)
p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) +
geom_boxplot(outlier.shape = NA, alpha = 0.1) +
geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) +
theme_ipsum_rc() +
facet_grid(trt~ .) +
scale_colour_ipsum() +
scale_fill_ipsum() +
theme(legend.title = element_blank()) +
labs(x = 'Exposure week', y = yvar)
p
Models with week
Olympia oyster
dat <- wtmods %>%
filter(species == 'Olympia') %>%
filter(var == yvar)
dat %>%
filter(species == 'Olympia') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Tissue weight (g) vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0 |
0 |
4 |
62.14 |
0.56 |
0.69 |
| week |
0 |
0 |
2 |
60.84 |
1.40 |
0.26 |
dat %>%
filter(species == 'Olympia') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Tissue weight (g) vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.00 |
0 |
4 |
52.47 |
0.60 |
0.66 |
| week |
0.00 |
0 |
2 |
53.14 |
1.36 |
0.27 |
| trt:week |
0.01 |
0 |
8 |
54.11 |
0.62 |
0.76 |
dat %>%
pull(plomod)
Pacific oyster
dat <- wtmods %>%
filter(species == 'Pacific') %>%
filter(var == yvar)
dat %>%
filter(species == 'Pacific') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Tissue weight (g) vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0 |
0 |
4 |
53.56 |
0.74 |
0.57 |
| week |
0 |
0 |
2 |
50.79 |
0.88 |
0.42 |
dat %>%
filter(species == 'Pacific') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Tissue weight (g) vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.00 |
0 |
4 |
42.28 |
0.92 |
0.46 |
| week |
0.00 |
0 |
2 |
41.75 |
0.96 |
0.39 |
| trt:week |
0.01 |
0 |
8 |
42.55 |
1.28 |
0.28 |
dat %>%
pull(plomod)
Whole weight
Raw data
yvar <- 'Whole weight (g)'
toplo <- allexp %>%
filter(var %in% yvar)
p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) +
geom_boxplot(outlier.shape = NA, alpha = 0.1) +
geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) +
theme_ipsum_rc() +
facet_grid(trt~ .) +
scale_colour_ipsum() +
scale_fill_ipsum() +
theme(legend.title = element_blank()) +
labs(x = 'Exposure week', y = yvar)
p
Models with week
Olympia oyster
dat <- wtmods %>%
filter(species == 'Olympia') %>%
filter(var == yvar)
dat %>%
filter(species == 'Olympia') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Whole weight (g) vs week + treatment
| |
Df |
Deviance |
Resid. Df |
Resid. Dev |
| NULL |
NA |
NA |
305 |
23.08 |
| trt |
4 |
0.61 |
301 |
22.47 |
| week |
2 |
1.05 |
299 |
21.42 |
dat %>%
filter(species == 'Olympia') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Whole weight (g) vs week x treatment
| |
Df |
Deviance |
Resid. Df |
Resid. Dev |
| NULL |
NA |
NA |
305 |
23.08 |
| trt |
4 |
0.61 |
301 |
22.47 |
| week |
2 |
1.05 |
299 |
21.42 |
| trt:week |
8 |
1.03 |
291 |
20.39 |
dat %>%
pull(plomod)
Pacific oyster
dat <- wtmods %>%
filter(species == 'Pacific') %>%
filter(var == yvar)
dat %>%
filter(species == 'Pacific') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Whole weight (g) vs week + treatment
| |
Df |
Deviance |
Resid. Df |
Resid. Dev |
| NULL |
NA |
NA |
290 |
7.68 |
| trt |
4 |
0.24 |
286 |
7.44 |
| week |
2 |
0.15 |
284 |
7.28 |
dat %>%
filter(species == 'Pacific') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Whole weight (g) vs week x treatment
| |
Df |
Deviance |
Resid. Df |
Resid. Dev |
| NULL |
NA |
NA |
290 |
7.68 |
| trt |
4 |
0.24 |
286 |
7.44 |
| week |
2 |
0.15 |
284 |
7.28 |
| trt:week |
8 |
1.07 |
276 |
6.22 |
dat %>%
pull(plomod)
Delta length
Raw data
yvar <- 'Delta length (cm)'
toplo <- allexp %>%
filter(var %in% yvar)
p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) +
geom_boxplot(outlier.shape = NA, alpha = 0.1) +
geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) +
theme_ipsum_rc() +
facet_grid(trt~ .) +
scale_colour_ipsum() +
scale_fill_ipsum() +
theme(legend.title = element_blank()) +
labs(x = 'Exposure week', y = yvar)
p
Models with week
Olympia oyster
dat <- wtmods %>%
filter(species == 'Olympia') %>%
filter(var == yvar)
dat %>%
filter(species == 'Olympia') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Delta length (cm) vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.07 |
0.02 |
4 |
45.09 |
2.28 |
0.08 |
| week |
0.05 |
0.05 |
1 |
45.08 |
5.73 |
0.02 |
dat %>%
filter(species == 'Olympia') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Delta length (cm) vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.07 |
0.02 |
4 |
40.94 |
2.15 |
0.09 |
| week |
0.04 |
0.04 |
1 |
40.95 |
5.41 |
0.03 |
| trt:week |
0.01 |
0.00 |
4 |
40.94 |
0.34 |
0.85 |
dat %>%
pull(plomod)
Pacific oyster
dat <- wtmods %>%
filter(species == 'Pacific') %>%
filter(var == yvar)
dat %>%
filter(species == 'Pacific') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Delta length (cm) vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.06 |
0.01 |
4 |
44.26 |
1.67 |
0.17 |
| week |
0.00 |
0.00 |
1 |
44.26 |
0.22 |
0.64 |
dat %>%
filter(species == 'Pacific') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Delta length (cm) vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.06 |
0.01 |
4 |
40.24 |
1.66 |
0.18 |
| week |
0.00 |
0.00 |
1 |
40.24 |
0.21 |
0.65 |
| trt:week |
0.03 |
0.01 |
4 |
40.24 |
0.90 |
0.47 |
dat %>%
pull(plomod)
Delta width
Raw data
yvar <- 'Delta width (cm)'
toplo <- allexp %>%
filter(var %in% yvar)
p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) +
geom_boxplot(outlier.shape = NA, alpha = 0.1) +
geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) +
theme_ipsum_rc() +
facet_grid(trt~ .) +
scale_colour_ipsum() +
scale_fill_ipsum() +
theme(legend.title = element_blank()) +
labs(x = 'Exposure week', y = yvar)
p
Models with week
Olympia oyster
dat <- wtmods %>%
filter(species == 'Olympia') %>%
filter(var == yvar)
dat %>%
filter(species == 'Olympia') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Delta width (cm) vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.06 |
0.02 |
4 |
44.07 |
1.36 |
0.26 |
| week |
0.04 |
0.04 |
1 |
44.07 |
3.36 |
0.07 |
dat %>%
filter(species == 'Olympia') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Delta width (cm) vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.06 |
0.01 |
4 |
40.07 |
1.25 |
0.30 |
| week |
0.04 |
0.04 |
1 |
40.07 |
3.08 |
0.09 |
| trt:week |
0.00 |
0.00 |
4 |
40.07 |
0.09 |
0.99 |
dat %>%
pull(plomod)
Pacific oyster
dat <- wtmods %>%
filter(species == 'Pacific') %>%
filter(var == yvar)
dat %>%
filter(species == 'Pacific') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Delta width (cm) vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.03 |
0.01 |
4 |
44.19 |
0.50 |
0.74 |
| week |
0.00 |
0.00 |
1 |
44.19 |
0.06 |
0.81 |
dat %>%
filter(species == 'Pacific') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Delta width (cm) vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.03 |
0.01 |
4 |
40.2 |
0.49 |
0.74 |
| week |
0.00 |
0.00 |
1 |
40.2 |
0.06 |
0.81 |
| trt:week |
0.05 |
0.01 |
4 |
40.2 |
0.89 |
0.48 |
dat %>%
pull(plomod)
Length rate
Raw data
yvar <- 'Length rate (cm/d)'
toplo <- allexp %>%
filter(var %in% yvar)
p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) +
geom_boxplot(outlier.shape = NA, alpha = 0.1) +
geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) +
theme_ipsum_rc() +
facet_grid(trt~ .) +
scale_colour_ipsum() +
scale_fill_ipsum() +
theme(legend.title = element_blank()) +
labs(x = 'Exposure week', y = yvar)
p
Models with week
Olympia oyster
dat <- wtmods %>%
filter(species == 'Olympia') %>%
filter(var == yvar)
dat %>%
filter(species == 'Olympia') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Length rate (cm/d) vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0 |
0 |
4 |
45.41 |
2.19 |
0.08 |
| week |
0 |
0 |
1 |
45.41 |
1.22 |
0.27 |
dat %>%
filter(species == 'Olympia') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Length rate (cm/d) vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0 |
0 |
4 |
41.53 |
2.16 |
0.09 |
| week |
0 |
0 |
1 |
41.54 |
1.20 |
0.28 |
| trt:week |
0 |
0 |
4 |
41.53 |
0.83 |
0.51 |
dat %>%
pull(plomod)
Pacific oyster
dat <- wtmods %>%
filter(species == 'Pacific') %>%
filter(var == yvar)
dat %>%
filter(species == 'Pacific') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Length rate (cm/d) vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0 |
0 |
4 |
44.22 |
1.85 |
0.14 |
| week |
0 |
0 |
1 |
44.22 |
4.38 |
0.04 |
dat %>%
filter(species == 'Pacific') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Length rate (cm/d) vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0 |
0 |
4 |
40.22 |
1.91 |
0.13 |
| week |
0 |
0 |
1 |
40.22 |
4.50 |
0.04 |
| trt:week |
0 |
0 |
4 |
40.22 |
1.30 |
0.29 |
dat %>%
pull(plomod)
Width rate
Raw data
yvar <- 'Width rate (cm/d)'
toplo <- allexp %>%
filter(var %in% yvar)
p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) +
geom_boxplot(outlier.shape = NA, alpha = 0.1) +
geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) +
theme_ipsum_rc() +
facet_grid(trt~ .) +
scale_colour_ipsum() +
scale_fill_ipsum() +
theme(legend.title = element_blank()) +
labs(x = 'Exposure week', y = yvar)
p
Models with week
Olympia oyster
dat <- wtmods %>%
filter(species == 'Olympia') %>%
filter(var == yvar)
dat %>%
filter(species == 'Olympia') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Width rate (cm/d) vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0 |
0 |
4 |
44.22 |
1.18 |
0.33 |
| week |
0 |
0 |
1 |
44.22 |
0.96 |
0.33 |
dat %>%
filter(species == 'Olympia') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Width rate (cm/d) vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0 |
0 |
4 |
40.21 |
1.10 |
0.37 |
| week |
0 |
0 |
1 |
40.21 |
0.90 |
0.35 |
| trt:week |
0 |
0 |
4 |
40.21 |
0.29 |
0.88 |
dat %>%
pull(plomod)
Pacific oyster
dat <- wtmods %>%
filter(species == 'Pacific') %>%
filter(var == yvar)
dat %>%
filter(species == 'Pacific') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Width rate (cm/d) vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0 |
0 |
4 |
44.21 |
0.30 |
0.88 |
| week |
0 |
0 |
1 |
44.22 |
10.89 |
0.00 |
dat %>%
filter(species == 'Pacific') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Width rate (cm/d) vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0 |
0 |
4 |
40.2 |
0.29 |
0.88 |
| week |
0 |
0 |
1 |
40.2 |
10.46 |
0.00 |
| trt:week |
0 |
0 |
4 |
40.2 |
0.56 |
0.69 |
dat %>%
pull(plomod)
Final length
Raw data
yvar <- 'Final length (cm)'
toplo <- allexp %>%
filter(var %in% yvar)
p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) +
geom_boxplot(outlier.shape = NA, alpha = 0.1) +
geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) +
theme_ipsum_rc() +
facet_grid(trt~ .) +
scale_colour_ipsum() +
scale_fill_ipsum() +
theme(legend.title = element_blank()) +
labs(x = 'Exposure week', y = yvar)
p
Models with week
Olympia oyster
dat <- wtmods %>%
filter(species == 'Olympia') %>%
filter(var == yvar)
dat %>%
filter(species == 'Olympia') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Final length (cm) vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.16 |
0.04 |
4 |
52.98 |
0.45 |
0.77 |
| week |
0.24 |
0.12 |
2 |
52.95 |
1.35 |
0.27 |
dat %>%
filter(species == 'Olympia') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Final length (cm) vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.20 |
0.05 |
4 |
44.90 |
0.54 |
0.71 |
| week |
0.26 |
0.13 |
2 |
44.91 |
1.44 |
0.25 |
| trt:week |
1.05 |
0.13 |
8 |
44.91 |
1.44 |
0.21 |
dat %>%
pull(plomod)
Pacific oyster
dat <- wtmods %>%
filter(species == 'Pacific') %>%
filter(var == yvar)
dat %>%
filter(species == 'Pacific') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Final length (cm) vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.20 |
0.05 |
4 |
52.96 |
0.56 |
0.70 |
| week |
0.06 |
0.03 |
2 |
52.91 |
0.33 |
0.72 |
dat %>%
filter(species == 'Pacific') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Final length (cm) vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.14 |
0.03 |
4 |
44.86 |
0.38 |
0.82 |
| week |
0.06 |
0.03 |
2 |
44.89 |
0.34 |
0.72 |
| trt:week |
0.73 |
0.09 |
8 |
44.89 |
1.02 |
0.43 |
dat %>%
pull(plomod)
Final width
Raw data
yvar <- 'Final width (cm)'
toplo <- allexp %>%
filter(var %in% yvar)
p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) +
geom_boxplot(outlier.shape = NA, alpha = 0.1) +
geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) +
theme_ipsum_rc() +
facet_grid(trt~ .) +
scale_colour_ipsum() +
scale_fill_ipsum() +
theme(legend.title = element_blank()) +
labs(x = 'Exposure week', y = yvar)
p
Models with week
Olympia oyster
dat <- wtmods %>%
filter(species == 'Olympia') %>%
filter(var == yvar)
dat %>%
filter(species == 'Olympia') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Final width (cm) vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.18 |
0.05 |
4 |
53.03 |
0.67 |
0.61 |
| week |
0.23 |
0.12 |
2 |
52.99 |
1.69 |
0.19 |
dat %>%
filter(species == 'Olympia') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Final width (cm) vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.17 |
0.04 |
4 |
44.97 |
0.61 |
0.66 |
| week |
0.23 |
0.12 |
2 |
44.99 |
1.72 |
0.19 |
| trt:week |
0.60 |
0.08 |
8 |
44.99 |
1.11 |
0.38 |
dat %>%
pull(plomod)
Pacific oyster
dat <- wtmods %>%
filter(species == 'Pacific') %>%
filter(var == yvar)
dat %>%
filter(species == 'Pacific') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Final width (cm) vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.06 |
0.02 |
4 |
53.00 |
0.31 |
0.87 |
| week |
0.11 |
0.06 |
2 |
52.93 |
1.07 |
0.35 |
dat %>%
filter(species == 'Pacific') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Final width (cm) vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
0.12 |
0.03 |
4 |
44.85 |
0.56 |
0.70 |
| week |
0.11 |
0.05 |
2 |
44.89 |
1.05 |
0.36 |
| trt:week |
0.38 |
0.05 |
8 |
44.89 |
0.91 |
0.52 |
dat %>%
pull(plomod)
Respiration
Raw data
yvar <- 'Respiration (umol/hr/g)'
toplo <- allexp %>%
filter(var %in% yvar)
p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) +
geom_boxplot(outlier.shape = NA, alpha = 0.1) +
geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) +
theme_ipsum_rc() +
facet_grid(trt~ .) +
scale_colour_ipsum() +
scale_fill_ipsum() +
theme(legend.title = element_blank()) +
labs(x = 'Exposure week', y = yvar)
p
Models with week
Olympia oyster
dat <- wtmods %>%
filter(species == 'Olympia') %>%
filter(var == yvar)
dat %>%
filter(species == 'Olympia') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Respiration (umol/hr/g) vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
16.68 |
4.17 |
4 |
51.04 |
0.34 |
0.85 |
| week |
9.67 |
4.83 |
2 |
55.99 |
0.40 |
0.67 |
dat %>%
filter(species == 'Olympia') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Respiration (umol/hr/g) vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
13.93 |
3.48 |
4 |
44.51 |
0.29 |
0.88 |
| week |
10.23 |
5.12 |
2 |
46.96 |
0.42 |
0.66 |
| trt:week |
64.95 |
8.12 |
8 |
47.14 |
0.67 |
0.72 |
dat %>%
pull(plomod)
Pacific oyster
dat <- wtmods %>%
filter(species == 'Pacific') %>%
filter(var == yvar)
dat %>%
filter(species == 'Pacific') %>%
pull(anomodnoint) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Respiration (umol/hr/g) vs week + treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
40.25 |
10.06 |
4 |
53.53 |
0.73 |
0.58 |
| week |
11.31 |
5.66 |
2 |
52.35 |
0.41 |
0.67 |
dat %>%
filter(species == 'Pacific') %>%
pull(anomod) %>%
.[[1]] %>%
knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
kable_styling(full_width = T, font_size = 14) %>%
HTML
Respiration (umol/hr/g) vs week x treatment
| |
Sum Sq |
Mean Sq |
NumDF |
DenDF |
F value |
Pr(>F) |
| trt |
37.54 |
9.39 |
4 |
45.50 |
0.69 |
0.61 |
| week |
9.84 |
4.92 |
2 |
44.45 |
0.36 |
0.70 |
| trt:week |
54.96 |
6.87 |
8 |
45.48 |
0.50 |
0.85 |
dat %>%
pull(plomod)
